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A  THEOREM  ON  TREND- FREE  BLOCK  DESIGNS 


Ching-Ming  Yen  Ralph  A.  Bradley 

University  of  flinnesota  at  Morris  Florida  State  University 

Morris,  Minnesota  Tallahassee,  Florida 


SUt*!ARY 


-J 

An  earlier,  referenced  report  considered  the  existence  of  trend- free 
block  designs  md  gave  some  general  results  for  general  block  designs  and 
some  special  results  for  corplete  block  designs.  This  report  presents  a 
new  existence  theorem  that  is  sufficiently  general  to  cover  broad  classes 
of  incomplete  block  designs  and  suggests  a  design  construction  method. 

Suppose  that  a  common  trend  affecting  treatment  responses  exists  over 
the  plots  within  blocks  of  a  block  design  with  the  usual  design  parameters, 
v,  b,  k,  r.  A  trend- free  block  design,  if  it  exists,  is  one  for  which 
treatments  are  assigned  to  plots  within  blocks  so  that  treatment  aad  block 
contrasts  are  orthogonal  to  specified  trend  components .  V.'e  prove  that , 
if  a  given  block  design  is  constructed  and  the  common  trend  is  linear,  the 
necessary  and  sufficient  condition  for  the  existence  of  a  trend- free  block 
design  with  the  same  incidence  matrix  is  that  %r(k  +  i)  be  an  integer. 
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Bradley  and  Yeh  (1930)  introduced  the  concept  of  trend- free  block  (TFB) 

designs  and  developed  a  necessary  and  sufficient  condition  for  their  existence. 

Yeh  and  Bradley  (1979)  established  alternative  existence  conditions,  in 

some  cases  for  special  trends  and  special  block  designs,  and  examined  some 

design  construction  procedures.  While  some  important  results  were  obtained 

for  complete  block  designs,  extensions  to  incomplete  block  designs  were 

limited.  The  purpose  of  this  report  is  to  present  a  new  theorem  tliat  is 

very  generally  applicable  to  incomplete  block  designs. 

In  review,  consider  an  experimental  situation  with  v  treatments 

applied  to  plots  arranged  in  b  blocks  of  size  k  s  v.  Each  plot  receives 

only  one  treatment  and  each  treatment  occurs  at  most  once  in  a  block. 

Plots  in  each  block  are  arranged  in  the  same  m- dimensional  array,  indexed 

by  vectors  of  positive  integers,  t  =  (tp  ...,  t^) ,  tu  *  1,  ...,  su, 

u  *  1,  ...,  n,  where  su  is  the  number  of  plot  positions  in  the  utJl 

dimension,  n  s  =  k.  A  polynomial  trend;  common  to  all  blocks ,  is 
u 

assumed  to  exist  over  the  plots  in  a  block  and  to  be  a  function  of  the 
plot  position  t. 

The  classical  model  for  general  block  designs  is  extended  through 
the  addition  of  trend  terns.  Die  trend  function  is  expressed  as  a  linear 
corrbination  of  n- dimensional,  orthogonal  polynomials  of  the  form. 


»>  rt)  =  „  .(v 


. t  )  =  n  4>  (t  ) , 
rr  a  u'  ' 

u=l  u 


where  ^  (tu)  is  a  one -dimensional ;  orthogonal  polynomial  of  degree  on 
the  integers,  1,  ...,  su<  The  model  is 


■  ■  *  ilM  * 6;  *  *  ■«- 


(2) 
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j  =  1,  ....  b,  t  =  (tr  ....  tm),  tu=  1,  ...,  su,  u=  1,  ...,  m,  where 
y.  is  the  observation  on  plot  position  t  of  block  j,  w,  t.  and  6-  are 

•''V  ^  J 

respectively  the  usual  mean,  treatment  and  block  parameters,  £  6  4>  (t)  is 

jjeA  ®  ® 

the  trend  effect  on  plot  t,  not  dependent  on  the  particular  block  j ,  with 
0q  being  the  regression  coefficient  of  $  (t)  and  A,  an  index  set  of  p, 
m- dimensional,  non-zero  vectors  of  tlie  form  a,  p  <  k,  and  the  e-  are 
random  errors  assumed  to  be  i.i.d.  with  zero  means.  Designation  of  the 
treatment  applied  to  plot  (j ,  t)  is  effected  through  indicator  variables , 

6 j t  =  1  or  0  as  treatment  i  is  or  is  not  on  plot  ( j ,  t) ,  i  =  1 ,  . . . ,  v. 

A  TFB  design  exists  when  it  is  possible  to  choose  the  so  that  the 
estimators  of  treatment  and  block  contrasts  arc  unaffected  by  the 
presence  or  absence  of  the  trend  terms  in  (2).  For  a  TFB  design,  the 
appropriate  treatment  and  block  sums  of  squares  in  analysis  of  variance 
have  the  same  algebraic  forms  in  the  presence  or  absence  of  the  trend 
terms. 

Let  A j ,  j  =  1 ,  . . . ,  b ,  be  the  k  *  v  matrix  with  6 j t  in  row  t  and  column 
i  and  let  <t>  be  the  k  x  p  matrix  with  <J>  (t)  in  row  t  and  column  a. 

Further,  let  XJ.  =  (Aj,  ...,  and  X0  =  1^  8  <j>,  where  1^  is  the 
b-dinensional  colum  vector  with  unit  elements  and  B  9  C  is  the  Kronecker 
product  of  B  and  C.  Bradley  and  Yell  (1980,  Theorem  3.1)  show  that  a  block 
design  is  trend- free  (to  tlie  trend  tenns  specified  in  the  model)  if  and 
only  if 

(3)  X'  X  =  0, 

v  J  ~T  ~6  ~ 

or  equivalently,  if  and  only  if 

(4)  a;  ♦  -  o, 

l) 

where  A^  =  J  &■. 

j-l"3 


In  this  report,  attention  is  limited  to  consideration  of  TFB  designs 
when  =  ...  =  r^  =  r,  k  s  v,  r,  k  >  1,  m  =  1.  These  designs  will  be 
designated  as  TF^B(v,  b,  k,  r)  designs,  designs  with  treatments  so 
arranged  within  blocks  that  they  are  trend- free  to  a  common,  one* dimensional 
trend  over  the  k  plot  positions  within  the  blocks.  When  $  =  9^.  the 
common  trend  is  linear;  when  9  =  (9. ,  $9)  ,  the  common  trend  is  quadratic. 

Theorem  3.3  of  Yen  and  Bradley  (1979)  is  the  key  to  the  new  theorem 
of  this  report  and  is  -estated: 


Theorem  1.  If  <j>  =  (9.. ,  . . .  .  4>  )  under  model  (2),  then  a  TF  B(v,  b.  k,  r) 

-  ~  ~p  9 


design  exists  if  and  only  if  there  exists  a  v  x  b  matrix  W  with  non-negative 
elements  such  tliat 

(i)  each  coluim  of  W  lias  the  integers  1,  ...»  k  as  elements  along 
with  (v  -  k)  zero  elements,  and, 


th 


(ii)  for  any  a  *  1,  ...  p,  the  sun  of  the  a  powers  of  elements 

k 

for  any  single  row  of  W  is  Sa(k»  r)  =  J  i“. 


i=l 


Hie  theorem  to  be  proved  is . 


Tneorem  2.  If  9  *  9,  and  a  connected  block  design  exists  with 

Tr“  »v  <vl 


parameters  v,  b,  k  and  r,  tlie  design  may  be  converted  into  a  TF  D(v,  b,  k,  r) 

21 

design  throu^i  exchanges  of  plot  positions  for  treatments  within  blocks  if 


and  only  if  5ir(k  +  1)  is  an  integer. 


Necessity  follows  at  once  from  Theorem.  1.  Proof  of  sufficiency  requires 
a  series  of  lemmas  given  below.  Theorem  2  irplies  Theoren  4.2  (Yell  and 
Bradley,  1379)  when  k  =  v  and  application  to  various  classes  of  incomplete 
block  designs  is  obvious. 
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The  key  to  proof  of  sufficiency  for  Theorem  2  is  the  matrix  of 
Theorem  1.  The  typical  element  w^  of  1/  is  the  plot  position  of  treatment 
i  in  block  j  if  f  0  and,  if  w.^  =  0,  treatment  i  does  not  occur  in 
block  j.  Each  colum  of  W  consists  of  the  integers  1,  . ..,  k  along  with 
v  -  k  zeros  and  each  row  of  W  contains  r  non- zero  elements.  The  first  and 
obvious  lemma  summarizes  the  situation. 


Lemma  1.  If  W  represents  any  block  design  with  parameters 
v,  b,  k  and  r  with  specified  plot  positions  for  treatments  within  each 
block,  then 

(i)  (The  minimum  row  sum  of  W)  5  %r(lc  +  1) ,  and 

(ii)  equality  holds  in  (i)  if  and  only  if  W  represents  a 

TF  B(v,  b,  k,  r)  design. 

Zl 

Lemma  1  provides  a  basis  for  the  desired  proof  of  sufficiency.  If 

exchanges  of  non- zero  elements  are  made  within  columns  of  W,  a  new  design 

with  the  designated  parameters  results  but  row  sums  in  V!  may  change.  If 

such  exchanges  result  in  equal  row  sums  of  ^r(k  +  1) ,  an  integer  by 

assumption  in  Theorem.  2,  then  a  TF  B(v,  b,  k,  r)  design  is  produced 

XI 

(Theorem  1) .  A  matrix  U  may  have  several  rows  with  minimum  row  sum. 

Definition  1.  A  matrix  VI,  and  its  corresponding  block  design,  con 
be  AjupKovzd.  if,  through  exchanges  of  non-zero  elements  in  colums  of  U, 
the  minimum  row  sura  of  the  resultant  matrix  is  not  decreased  from  that  of 
V!  and  the  nurber  of  rows  with  the  minimum  row  sum  of  W  is  reduced. 


Remarks  1.  If  4r(k  +  1)  is  an  integer  and  if  VI,  representing  a  block 
design  of  the  class  under  consideration,  is  such  that  it  can  be  improved 


through  a  sequence  of  in^roved  matrices 


\I  ,  VI  having  minimum 
~n’  ~n 
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row  sura,  hr(lc  +  1) ,  then  n  is  finite  and  V/  represents  and  specifies  a 

~n 

TF.  B(v,  b,  k,  r)  design. 

xi 


Some  exchange  notation  is  needed.  Wc  write 

(i)  w.  .  -*■  w. . .  to  mean  that  w. .  +  1  =  w. . . ,  i  f  i' ,  v. , .  *  0, 

ii  13  13  13  ’i*3 

(ii)  5  vc  to  mean  that  (i)  holds  and  the  exchange  is  made 
yielding  a  new  !7  matrix,  and 

(iii)  w.  .  5  w.  .  |w.  .  §  w.  .  to  mean  tliat  the  first  exc’n 

11J1  12J1  12J2  SJ2 

made  in  IV  followed  by  the  second  exchange  in  the  matrix  resulting  from  the 

first,  ^  f  j2. 

More  generally,  we  write 


excnange  is 


(5)  w.  .  Sw.  .  |w  ,  5  w.  -  | 

Vl  12J1  l2J2  VI 


lwi  a  +  w.  i 
^•’n  ^tt+l^n 


to  mean  a  sequence  of  exchanges,  each  based  on  the  h  matrix  resulting 
from  the  previous  exchange,  jt  t  j  p  t  =  1,  n. 

Lerma  2.  Given  the  sequence  of  exchanges  (5)  with  n  >  1,  the 
resultant  IV  matrix  is  such  that, 

(i)  If  i^  f  i  p  the  i^-tli  row  sum  is  increased  by  1  (relative  to 
the  original  matrix  U),  the  in+^-th  row  sum  is  decreased  by  1,  and  all  of 
the  remaining  row  sums  are  unchanged,  or 

(ii)  If  i^  =  in+1,  all  row  sums  are  unchanged. 

Definition  2.  If  the  sequence  of  exchanges  (5)  is  possible  with 

i.  /  i  ,  we  say  that  row  i,  can  fizadi  row  i  in  the  initial  matrix  VV. 

1  n+1  i  n+1  ~ 

If  A  and  3  are  two  matrices  with  rows  consisting  of  two,  non-empty,  disjoint 
subsets  of  the  rows  of  W.  we  say  that  A  can  fimdi  B  if  a  row  in  A  can  reach 
a  row  in  B. 

! 
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Lenna  3.  If  A  and  B  contain  all  row  suns  in  W  and  if  W  represents 
a  connected  design  wit}i  parameters  v,  b,  k  and  r,  then  A  can  reach  B 
through  one  exchange  if  the  average  row  sum  of  A  does  not  exceed  hr(k  +  1). 


Proof.  The  proof  is  by  contradiction.  Suppose  that  A  cannot  reach  B 
through  one  exchange. 

Let  n.  be  the  number  of  non- zero  elements  in  the  j-th  column  of  A, 

J 

j  =  1,  ....  b.  If  n^  f  0,  the  non- zero  elements  in  column  j  of  A  must  be 
k,  (k  -  1),  ...,  (k  -  n j  +  1) >  the  largest  n^  numbers  in  the  set  {1,  ...»  k}. 
If  the  number  of  rows  in  A  is  £  k,  0  <  n^  <  k  for  at  least  one  j  since  the 
design  is  connected.  It  follows  that  the  average  row  sum  in  A  is 


b  b  _  b  7  b 

(6)  7  nj(21i  "  nj  +  nj  =  r(*  +  '  7^  Xl  V X1  "j-1 

b  b 

>  r(k  +  h)  -  y(  l  n.k/  £  n  )  =  I(k  +  1) , 
1  j=l  J  j=l  J 


the  inequality  in  (6)  being  strict  because  the  design  is  connected.  The 
condition  of  the  lemma  is  violated  and  the  proof  is  complete. 

Lemma  4  follows  from  Lemma  2  and  Lemma  5  presents  a  condition  for 
continuation  of  the  exchange  process: 


Lemma  4.  If  rows  i  and  i'  in  V/  have  respectively  the  minimum  row  sum 
n  and  a  row  sum  s  m  +  2,  W  can  be  improved  if  row  i  can  reach  row  i! . 


Leti.ia  5.  If  w. 


Vi 


w.  .  is  possible  in  U  and  if  the  resultant  matrix 
12^1 


does  not  represent  a  TF  B  design  but  does  have  row  i?  with  minimum  row 

xi  1 

sum,  there  exists  a  row  i^  and  a  column  j 2 »  *3  t  *2 '  -*2  ^  ^1’  sucl1  tnat 


v/ 


ij) 


w. 


2J  2 


13?2 


Proof.  Assume  that  the  clained  result  is  false.  Let  W*  with 

F 

typical  element  w'v.  result  after  w.  .  w.  .  in  W.  Then  row  of 

13  11J1  Vl  ~ 

W*  cannot  reach  any  other  row  through  a  column  different  from  j^.  This 

means  that  w*  .  is  now  at  least  1  and  all  of  the  remaining  (r  -  1) 

12J1 

positive  elements  in  row  i.  of  V/*  must  be  k.  Therefore,  the  i„-th  row 

L.  ~  ^ 

sun  of  W*  s  (r  -  l)k  +  1  s  Jjr(k  +  1)  since  r,  k  >  1.  A  contradiction  is 
established  by  Lenina  1  and  the  proof  follows. 

Notice  that,  if  i^  =  i^  in  Lemma  5,  each  row  sum  in  the  matrix  resulting 
E  E 

from  w.  ■  +w.  -I  w.  .  +  w.  •  matches  the  corresponding  row  sum  of  W 

1131  Vl  1232  15J2 

by  Lenina  2. 


Proof  of  Sufficiency,  Theorem  2.  If  tlie  matrix  W  corresponding  to 

the  block  design  of  Theorem  2  does  not  represent  a  TF  B  design,  all  rows  of 

Zl 

W  may  be  collected  into  three  matrices,  Rfl,  Rp  and  such  that  R() 

consists  of  rows  with  minimum  sum  n,  of  rows  with  sums  (m  +  1),  and 

of  rows  witli  suns  s  (m  +  2).  We  show  that  Rq  can  alv;ays  reach 

Therefore,  by  Lenina  4,  V/  can  be  improved  and  a  TF  B(v,  b,  k,  r)  design  can 

£l 

be  obtained  by  Remark  1. 

It  is  easy  to  show  that,  when  the  row  sums  of  U  are  not  all  equal,  R^ 
may  or  nay  not  exist,  while  Rq  and  always  exist.  V.hen  R^  does  not  exist, 
we  may  set  A  =  Rq  and  B  =  and  use  Lemma  3  to  sliow  that  Rq  can  reach  I^. 

When  R-.  exists,  we  sliow  that  R,,  can  reach  R0  bv  contradiction. 


Assume  that  Rq  cannot  reach  It,.  Set  A  =  Rg  and  B  =  and  use 

~2 

Lemma  3  again,  since  Rq  cannot  reach  >  Rq  can  reach  R^  througli  one 

£ 

exchange,  say  w.  .  -*■  w.  .  ,  yielding  a  new  matrix  \L .  Row  i~  in  VL 

Vl  12J1  ~ 
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has  sun  m,  the  minimun  row  sum  in  ,  and,  by  Lerma  5,  a  further  exchange, 

E 

w  •  •  w,  i  /  j  9  t  j  i  >  nay  be  made  and  yields  a  new  matrix  W7 .  After 
2J2  i3-)2  L  1  ~L 

the  two  exchanges,  row  i^  of  has  sim  m  or  (m  -  1)  depending  on  whether  row 

ij  was  in  R^  or  R^,  in  either  case,  row  has  minimum  row  sum  in  iv^* 

Lemma  5  is  applied  again  and  an  additional  exchange  may  be  made. 

A  successive  exchange  process  has  been  established  that  continues 

without  ending  if  cannot  reach  •  But  the  number  of  possible  exchanges 

and  orders  of  exchanges  in  columns  of  W  is  finite  and  a  closed  loop  of 

exchanges  like  (5)  with  i^+^  =  i.  must  result.  Rows  of  U  in  the  loop  muse 
~ ~0 

come  from  ^  because  otherwise  fL  would  reach  JL. .  The  loop  is  closed  in 
~1 

the  sense  that  rows  in  tlie  loop  cannot  reach  t!ie  remaining  rows,  if  any,  in 

V 

^  and  cannot  reach  rows  in  P^,  again  because  then  RQ  would  reach 
~1 


A  contradiction  lias  resulted  and  is  established  through  use  of  Lemma  3. 


Let  A  be  the  subset  of  rows  of 


in  tlie  closed  loon  and  let  B  be  the 


remaining  rows  of  W.  It  is  clear  that  the  average  row  sum  in  A  does  not 
exceed  br(k  +  1)  and  A  should  be  able  to  reacli  B.  Theorem  2  has  been 
established. 


Remark  2.  Theorem  2  has  irplications  for  rank-order  statistics.  If 
treatments  are  ranked  within  blocks  of  the  connected  block  design  and 
treatment  rank  sums  are  considered  instead  of  position  suns,  the  theorem 
establishes  conditions  under  which  the  treatments  may  have  equal  rank  sums. 

Let  a  balanced  incomplete  block  (BIB)  design  free  of  a  trend  described 
by  (j)  be  designated  as  a  TFiBIB(v,  b,  k,  r.  A)  design,  •  ‘any  TF,  BIB(v,  b,  k,  r, 

~  9  £i 

designs  exist,  Raghavarao  (1971,  Table  5.10.1)  gives  solutions  for  BIB 
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designs  with  v  b  s  100,  r,  k  <  15.  dinety-one  designs  are  listed  of  which 
only  fifteen  do  not  yield  TF  BID  designs,  namely  numbers  1  6,  14,  IS,  24, 

♦l 

27,  35,  36,  46,  43,  52,  53,  69,  72,  and  S3. 

Symmetric  balanced  incorplete  block  designs,  for  which  v  =  b  and  hence 
k  =  r,  free  of  a  trend  described  by  may  be  designated  as  TF  SBIB(v,  k,  x) 
designs,  where  X  has  the  usual  meaning.  It  is  obvious  that  SBIB  designs  yield 
TF  SBIB(v,  k,  x)  designs  by  Theorem  2.  A  Youden  Square  design  is  a  special 

<vl 

case  of  an  SBIB  design  in  a  rectangular  array  with  the  columns  forming  a 
C3  design  and  the  rows  an  SBIB  design.  Youden  Squares  are  the  incomplete 
block  design  counterpart  to  Latin  Squares  in  complete  block  designs. 

Theorems  on  trends  in  one  dimension  of  higher  order  than  the  first 
are  difficult  to  develop.  An  example  of  a  TF^BIBO,  b,  k,  r,  x)  Awign  is  given 
by  Yeh  and  Bradley  (1979 j  when  $  =  (<j>. ,  <*>,,) ,  a  quadratic  trend,  and 

~  'n.J. 

v  *  5,  b  =  10,  k  =  3,  r  »  6,  x  =  3. 

The  method  of  proof  of  Theorem  2  suggests  a  possible  exchange  method 
for  derivation  of  TF^  B(v,  b,  k,  r)  designs  when  they  exist  that  should  be 
feasible  on  a  computer. 
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